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1 INTRODUCTION 



The recent advent of large-scale galaxy surveys have revolu- 
tionized the field of observational cosmology. Deep spectro- 
scopic surveys allows us to witness the young Universe when 
the building blocks of present-day galaxies are forming, and 
some of the well-known properties of the local galaxy popu- 
lations are about to realize. The enormous amount of data 
gathered by wide area surveys produce galaxy samples with 
exquisite statistical precision, which makes it possible to sin- 
gle out the most fundamental properties that govern the 
physics of galaxy formation from the medley of observables. 

Equally impressive has been the progress in the the- 
oretical understanding of the structure formation in the 
Universe. Techniques such as direct numerical simulations 
and semi-analytic models can now reproduce the observed 
properties of galaxies, such as the luminosity function 
and 2-point correlation function, color, mass-to-light ra- 
tios over large ranges of e nvironments and cosmic epochs 
ilKauffmann et al.lll999allbl ; ISpringel et"aT]|200ll ; ICole et all 
l200d ; IWhite fc ReesHl978f ). 

Yet another approach, the so-called halo model, which 
is phenomenological in nature, has enjoyed popularity over 
the recent years. An essential ingredient of this method is 
the halo occupation distribution (HOD), which refers to the 
way galaxies (or substructures of dark matter halos) "pop- 
ulate" dark matter halos. In general, an HOD description 
includes the mean number of galaxies per halo A'' as a func- 



ABSTRACT 

We investigate statistical properties of LRGs in a sample of X-ray selected galaxy 
clusters at intermediate redshift (0.2 < z < 0.6). The LRGs are selected based on 
carefully designed color criteria, and the cluster membership is assessed via photo- 
metric redshifts. As clusters and LRGs are both viewed as promising tracer of the 
underlying dark matter distribution, understanding the distribution of LRGs within 
clusters is an important issue. Our main findings include: 

1. The halo occupation distribution of LRGs inside our cluster sample is N(M) = 
k x (M/10 14 ) Q where a = 0.620 ± 0.105 and k = 1.425 ± 0.285 assuming a Poisson 
distribution for N(M). 

2. The halo occupation distribution of LRGs (N(M)) and the satellite distribution 
of LRGs (N — 1(M)) are both consistent with being Poisson. To be more quantitative, 
we find Var(N)/(N) = 1.428 ± 0.351 and Var(N - 1)/(N - 1} = 1.823 ± 0.496. 

3. The radial profile of LRGs within clusters when fitted with a NFW profile gives 
a concentration of 17.5^4 3 (6.0^'g) including (excluding) BLRGs (Brightest LRGs). 

We also discuss the implications of these observations on the evolution of massive 
galaxies in clusters. 

tion of halo mass, the probability distribution that a halo 
of mass M contains N galaxies P(N\M), and the relative 
distribution (both sp atial and velocity) of galax ies and dark 
matter within halos (|Berlind fc Weinb erg 2002). 



The halo model formalism allows fast exploration of a 
wide range of HODs; an HOD that reproduces the observed 
clustering properties and luminosity function of galaxies 
can be further studied to reveal the physical processes that 
lead to galaxy formation and understanding of cosmolog- 
ical parameters. Examples of using halo model formalism 
to reproduce observables in order to reveal parameters in 
cosmology, galaxy ev olution and formation includes (e.g . 



Abazaiian et alj 120051; IWhite et al.1 120071; lYoo et alj l20od ; 
Zheng fe Weinberdl2007l ; iKulkarni et al.ll2007h ~ 



Despite the success in both observational and theoret- 
ical sides, there remains some unsolved problems regarding 
the formation of the massive, (usually) early type, galaxies. 
These galaxies appear " red and dead" , with the majority of 
the stars forming at high redshift (z >~ 2) and evolving pas- 
sively since. Within the cold dark matter (CDM) paradigm, 
in which massive galaxies are built by smaller galaxies via 
mergers in the late times, mergers between gas poor systems 
( "dry" mergers) seem to be a promising route to form giant 
galaxies. Observationally, however, the overall importance 
of dry mergers is still under heated debate. 

Luminous Red Galaxies (LRGs) are massive galaxies 
composed mainly of old stars, with little or no on-going star 
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formation. They demonstrate very consistent spectral en- 
ergy distribution (SED). Their SEDs mainly consist of old 
star spectrum, most notably for the 4000A break. This al- 
lows one to pho tometrically determine the ir redshifts fairly 
accurately (see IPadmanabhan et alj|2005h . With the accu- 
rate photometric redshifts of LRGs, one can probe a larger 
volume of the universe, thus giving better constraints on the 
formation of massive galaxies. By studying the HOD of the 
LRGs, we aim to provide a simple quantitative description 
of these galaxies in massive dark matter halos, which will en- 
able direct comparison with predictions of galaxy formation 
models. 

Here we aim to provide observational constraints 
on the HOD of the LRGs based on a sample of 47 
inter mediate-redshift clu sters from the ROSAT 400d sur- 
vey l|Burenin et al-l feoOG). with photometric data from the 
Sloan Digital Sky Survey fSDSSi lstoughton et al1l2002f ). Us- 
ing X-ray properties of these clusters to define the cluster 
center and estimate the cluster binding mass, we determine 
the mean halo occupation number N as a function of mass 
from ~ 1 x 10 14 M Q to ~ 8 x 10 14 M© and also investigate 
the LRG distribution and luminosity distribution within the 
clusters. 

In fj2] we briefly describe the X-ray cluster catalog that 
we utilize and the construction of SDSS LRG sample. In 
Sj3] we present our method and findings on the LRG distri- 
butions within the clusters and the mean halo occupation 
number. We discuss what is a good mass tracers and evolu- 
tion of massive galaxies in ^5] Possible systematics that may 
affect our results are discussed in §3] 

Throughout the paper we assume the cosmological pa - 
rameters to be the WMAP values l|Spergel et alj [2006): 
f2 m /i 2 = 0.1277, h = 0.732 and the Hubble parameter 
H = 73/i73 km s _1 Mpc" 1 . 



2 DATA 

2.1 Cluster Sample 

Our cluster sample is drawn from th e 400 square degree 
ROSAT PSPC Galaxy Cluster Survey jBurenin et alj|2006h 
(hereafter the 400d su rvey) , which is an ext ension of the 160 
square degree survey l|Vikhlinin et al.lll998l ). The survey de- 
tects extended X-ray sources in archival ROSAT PSPC im- 
ages down to a flux limit of 1.4 x 10 -13 ergs _1 cm -2 , with 
extensive optical spectroscopic follow up. Out of the 266 
clusters detected in the survey, 47 lie within the redshift 
range 0.2 < z < 0.6 and are covered by SDSS DR5. The 
redshift range is chosen to be consistent with the photo- 
metric cuts designed to select a homogeneous LRG sample 
across a wide range in cosmic epochs. 

The cluster catalog from the 400d survey provides esti- 
mates of cluster center, redshift, and X-ray luminosity Lx, 
which is used to estimate the cluster mass. Some of the basic 
information of the clusters in our sample is given in Table 1. 



2.2 LRG Data from Sloan Digital Sky Survey 

The Sloan Digital Sky Survey has taken ugriz CCD 
images of 10 4 deg 2 of the high-latitude sky. A ded- 
icated 2.5m telescope at Apache Point Observatory 



images the sky in 5 bands between 3000A and 10000A 
l|Fukugita et al.l Il996l) usin g a dr ift-scanning, mosai c 
CCD camera jGunn et all 1 19981 ; iGunn et alj I 20061 ). 
detecting objects to a flux limit of r ~ 22.5 mag. 
The survey selects 10 6 targets for spectros copy, most 
of t he m galaxies with r < 17.77 mag (|Gunn et all 



19981; York. Adelman. Anderson. Anderson, et al. 2000l : 



Stoughton et alj " 20021 ). This spectroscopic follow-up uses 
two digital spectrographs on the same telescope as the 
imaging camera. Details of the galaxy sur vey are described 
in the galaxy target sel ection papers (jEisenstein et all 
l200ll ; IStrauss et alj |2002| ); other aspects of the survey 
are mainly described in the Early Data Release paper 
l|Stoughton et alj 2002 ). All the data processing, including 
astromet rv llPier et al. 120031) . source identificati on and pho- 
tomet ry dLupton et al " 1 1 1 
I2OO4J) , calibration 



200ll;lHogg et al.ll 200 it llvezic et al 



Fukugita et alj 
2002), spectroscopic target selection 



1996: 



200 it IStrauss et alj I2OO2I ; iRichards et~ai1 



troscopic fiber placement (IBlanton et al 



Smith et al 
(lEisenstein et al' 



2002), and spec- 



2003) are don e 



automatically via SDSS software (| Tucker, et alj 120061 ) 
The SDSS is well-underway, and has had six major releases 
l|Adelman-McCarthv et al.ll2007i ). 

We utilize the pho tometric LRGs from SDSS con- 
structed as described in IPadmanabhan et alj l|2005l . here- 
after P05). The LRGs have been very useful as a cosmologi- 
cal probe since they are typically the most luminous galaxies 
in the universe, thus they probe a larger volume than most 
other tracers. On top of this, they also have very regular 
spectral energy distributions and a prominent 4000A break, 
making photometric redshift estimation much easier than 
the other galaxies. We plot the color magnitude diagram for 
one of the cluster and show that the LRGs in the cluster 
are the bright red galaxies that follow nicely along the red 
sequence (see Fig[T}. 

Our selection criteria are based on the spectros copic 
selection of LRGs described in lEisenstein et all (|200ll ). ex- 
tended to lower apparent luminosities (P05). We select 
LRGs by choosing galaxies that both have colors consistent 
with an old stellar population, as well as absolute luminosi- 
ties greater than a chosen threshold. The first criterion is 
simple to implement since the uniform SEDs of LRGs im- 
ply that they lie on an extremely tight locus in the space 
of galaxy colors; we simply select all galaxies that lie close 
to that locus. More specifically, we can define three (not 
independent) colors that describe this locus, 



c_l = (r-i) -0.25(3 -r) -0.18 , 

d ± = (r-i)- 0.125( ff - r) , 

C || = 0.7(3 - r) + 1.2(r - i - 0.18) , 



(1) 



where g, r, and i are the SDSS model magnitudes in these 
bands respectively. We now make the following color selec- 
tions, 



Cut I: I cj_ |< 0.2; 

Cut II : d ± > 0.55, g - r > 1.4, 

as well as the magnitude cuts 

Cut I: rp c tro < 13.6 + C| |/0.3; 

rpetro < 19.7; 
Cut II : i<18.3 + 2d ± , 



(2) 
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Figure 1. Color- magnitude diagram of objects toward the field 
of Cluster 142. Those satisfying cuts I and II (see Eqn. [5]and[3]l 
are shown as crosses and squares respectively. The red circles are 
selected LRGs and the blue dots are objects detected in SDSS 
photometrically. As shown, the selected LRGs lie very systemat- 
ically along the red sequence of the cluster. 



i < 20. 



(3) 



Making two cuts (Cut I and Cut II) is convenient since 
the LRG color locus changes direction sharply as the 4000A 
break redshifts from the g to the r band; this division divides 
the sample into low redshift (Cut I, z < 0.4) and high red- 
shift (Cut II, z > 0.4) samples. More details of these color 
selection criteria are thoroughly described in P05. 

We do however apply slightly different cuts than those 
adopted in P05: we limit our samples to sky regions where 



E(B-V) < 0.08 

and data taken under seeing condition of 
FWHM < 2.0". 



(4) 



(5) 



These cuts in extinction and seeing are applied simply by 
excluding areas at which the galaxy overdensity drops sig- 
nificantly. Furthermore, there are a few regions in SDSS that 
have 60% more red objects and less blue objects; we decide 
to throw away these regions. 

We slice our LRG sample into two redshift bins: 0.2 < 
Zphoto < 0.4 and 0.4 < z v hoto < 0.6. We also regularized 
our redshift distribution as described in P05. For our sam- 
ple, we have 855534 galaxies, covering 2,025,731 resolution 
10 HEALpix pixels, each with area of 11.8 arcmin 2 , giving 
0.422334 gal/pix. 

We then estimate the photometric redshift of these 
LRGs with the algorithm developed by P05. The typical 
uncertainty of the photo-z's is 8 Z = er z /(l + z) » 0.03 (see 
P05). 



3 ANALYSIS 
3.1 Method 

We estimate the cluster virial mass M200 = (47r/3)r| 00 x 
200p c from the X -ray luminosity using the ma ss-luminosity 
relation given bv lReiprich fc Bohringeil (|2002h 



log 



1A 2 L X (0.1 - 2.4 keV) 
/if 3 2 10 40 erg s- 1 



A + a log 



1.4M 20( 



(6) 



where A = —20.055 and a = 1.652. The radius r2oo is de- 
fined such that the enclosed mean overdensity is 200 times 
the critical density p c . The corresponding angular extent is 
0200- The mass-luminosity sca ling relation provides a mass 
estimate accurate to < 50% (jReiprich fc Bohringerl l2002h 
and a virial radius r2oo estimate accurate to 15%. 

As we now have the redshifts and positions of these 
clusters, we locate the LRGs as described in £|2.2l in each of 
these clusters. We look for LRGs that are within a cylinder of 
radius #200 and length of A z = 0.06 from the cluster center 
in both position and redshift space (i.e. zlrg = z c ± 0.03). 
We choose 8 Z — 0.03 since that is the typical la error on 
the LRG photometric redshift (P05) and z c is the cluster 
redshift. More discussion on the choice of cluster radius and 
8 Z will be described in 

Since we are relying on the photometric redshifts of the 
LRGs to find out whether a LRG sits in certain cluster or 
not, we take into account the effects of the following mecha- 
nisms that may lead to over-(or underestimate of the num- 
ber of LRGs in each cluster: 

i. LRG identific ation failure: There is an i dentification 
failure rate of ~ 1% (Padmanabhan et al. 2006|). This is the 
rate of which a LRG (photometrically chosen) is actually a 
star or a quasar after we get the spectra of the object. 

ii. Interlopers: There is a finite probability of finding 
LRGs inside the cluster purely by chance, we call these in- 
terlopers. We access the expected number of interlopers in 
each cluster by looking at the average number of LRGs in 
sky (2D projected) in the solid angle of radius = 6*200 of the 
cluster and the average probability of finding a LRG in red- 
shift range of z c ±S z where S z = 0.03 (as defined above). We 
can write down the expected number of interlopers ({Ni n t}) 
as: 



(N int ) = mv6 



P(z p )dzp 



(7) 



where P(z p ) is the normalized (photometric) redshift distri- 
bution of LRGs, n is the 2D average LRG density. 

iii. Missing galaxies due to errors in photometric red- 
shift: As an LRG can be scattered out of the cluster (due to 
photoz error), we need to account for this process by look- 
ing at the probability of LRG having been photometrically 
determined to be outside of the cluster, but in fact has spec- 
troscopic redshift that falls within the range of the cluster: 



P(W- 



> 8 Z , \\z s —z c \ < 8 Z ) 



(F(z s )+B(z s ))dz s ,(8) 



z-min = max(0.05, z c — 5 Ze ,in) 
Zmax = mm(0.7, Zc + 8 Zs>in ) 



F(z) 



P(5,z)d5 
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Figure 2. Spatial distribution of LRGs in Cluster 142. The blue 
points represent all objects detected in the SDSS photometric 
survey. Those satisfying cuts I and II are shown as crosses and 
squares, respectively. The LRGs that have photo- 2 consistent with 
the cluster redshift (z = 0.353) are represented as circles. The 
black cross denotes the centroid position of the ICM. The green 
circle is the region encircled by the virial radius of the cluster. 

r+oo 

B{z) = / P(S,z)dS , (9) 

where P(8, z) is the probability of findi ng S (= z s — z p ) at 
z s , given by iPadmanabhan et al.l (|2005l ) and these are only 
characterized within the spectroscopic redshift range from 
2 = 0.05 to z — 0.7. 5z c ,in is the redshift range we allow a 
LRG to be a cluster member when we have its spectroscopic 
redshift, and this is set to be 0.01. 

We then calculated the corrected LRG counts in each 
cluster via the following: 

(N corr ) = {{N obs ) - (N mt ))/f(z p , z c , z.) , (10) 
f(z p ,z c ,z 3 ) = (l-P(\z P -z c \ > 8 z \\z s -z c \ < 6 x ))x(l+F)) ,(11) 

and F is the LRG identification failure rate. 

We list these corrected LRG counts in Table 1. 

To convert the observed magnitudes of the LRGs into 
the rest-frame luminosity at z — 0, we follow the evolution 
of a simple stellar population formed in a burst at z = 5, 
with solar metall icity and Salpeter in i tial m ass function, us- 
ing the model of iBruzual fc Chariot! (|2003l ). The LRGs are 
selected so that their present-day magnitude lies in the range 
—23.5 < M g < — 21 (roughly corresponding to 1-7L», where 
L» is the characteristic luminosity). 

For each cluster, we visually inspect the spatial and 
color distributions of LRGs with respect to all objects de- 
tected by SDSS. An example is shown for cluster 142. Per- 
haps not surprisingly, the spatial distribution of the LRGs 
seems concentrated towards cluster center (Fig. 0. 

A general scenario that has been painted about LRGs 
and clusters is that there is a massive red galaxy sitting right 
in the middle of the cluster. Then some other process may 



Figure 3. The distribution of LRGs in the clusters. The number 
of LRGs in each bin are normalized by dividing the number of 
LRGs in each bin by the total number of LRGs in all bins. 

sometimes bring in other massive red galaxies, but they will 
probably sink into the center over several dynamical times. 

Here, we actually have the ability to see if this scenario 
is true: we have the number of LRGs inside each of these 
clusters and we know where they are. Below we present re- 
sults on the spatial distribution of LRGs in clusters ( i)3.2fl 
and the halo occupation number ( i]3.3[l . and the LRG mul- 
tiplicity function ( >i3.5|l . 

3.2 Spatial Distribution of LRGs within Clusters 

We show the spatial distribution of LRGs within the clus- 
ters in Fig. B] Pre vious studies (e.g. I Jones fc Forman|[l983 : 
iLin fc Mohrll2004f ) have shown that brightest galaxies tend 
to lie at the center of the clusters. Here we test if this is true 
for the LRGs. We plot the distribution of brightest LRGs in 
each of the cluster alongside with their companions in each 
of the cluster (see Fig. [?}. One realizes that most (~ 80%) 
of the brightest LRGs resides within the inner 20% of the 
scaled radius of the cluster. Therefore, we are consistent with 
the picture of having brightest galaxies lying at the centers 
of clusters. However, there is a significant fraction of clusters 
that does not follow this rule. 

The question of whether the centers of intracluster gas 
coincide with the central LRGs (defined as the LRG closest 
to the centroid of the X-ray emitting gas) is also very impor- 
tant to the understanding of the formation of galaxies. We 
investigated the distribution of the central LRGs inside the 
cluster (see Fig. [4}. There are ~ 20% of the "central" LRGs 
which are not central at all. This may suggest a few scenar- 
ios, one being that the cluster is not relaxed enough for the 
central LRG to sit at the center of the gravitational potential 
(which is supposedly traced by the intracluster light). The 
centroiding of the clusters in X-ray is called into question, 
and we will address this in i|4] 

The profile of galaxies in cluster is a key ingredient to 
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Figure 4. (Above) The distribution of brightest LRGs (BLRGs) 
and the non-brightest LRGs in the clusters. As shown above, the 
BLRGs tend to lie at centers of the clusters, while the ones that 
are not BLRGs have a shallower radial distribution. (Below) The 
distribution of most central LRGs and the non-central LRGs in 
the clusters. The distribution of the central LRGs are very similar 
to the BLRGs. 



Figure 5. The distribution of LRGs in the clusters with the fit to 
the NFW profile. Blue line: we fit the surface density of the LRGs 
(including BLRGs) to a NFW profile and get a concentration of 
17.5^4 3 w ith X 2 = 4.29. Black line: we fit the surface density 
of the LRGs (excluding BLRGs) to a NFW profile and get a 
concentration of 6.0^'g with x 2 = 6.6. 



the halo model formalism. One would like to understand 
how statistically LRGs populate the clusters they are resid- 
ing. We try to fit the NFW profile here to the LRG sur- 
face density of stacked clusters in our sample and find that 
the concentration of the surface density to be I7.5I4 3 with 
X 2 = 4.29. We also plot the fitted profile in Fig. [5] We also fit 
the NFW profile to LRG surface density of stacked clusters 
without the BLRG (Brightest LRG in the cluster), and this 
gives a concentration of 6.0^'g with x 2 = 6.6. Both profiles 
have very simila r concentrati o n as t he K-band cluster pro- 
file discussed in iLin fc Mohrl (|2007l ) (see Figure [6}. Errors 
in r2oo determination do n ot affect the fit in a ny significant 
fashion as demonstrated in lLin fc Mohij (|2007l ) appendix . 

3.3 Halo Occupation Number 

As the halo occupation number of the clusters is a key ingre- 
dient to the halo model formalism, we investigate the num- 
ber of LRGs in these clusters as a function of their masses. 
As the size of our sample is not large and the mass estimate 
o the clusters are accurate to 30% — 50% only, one will have 
to be extra cautious in finding a fit for the average number 
of LRGs in the mass range of these clusters. We take the 
following approach, assuming two different models: 

N(M t ) = a x M t + k 
N(M t ) = kx Mt 

(12) 

where Mt is true value of cluster virial mass in 
a Poisson distribution of N(M t ) and two distributions for 
the probability finding Mt given Mi where Mi is the mea- 
sured mass of the i-th cluster (in same units as in M t ) , one 



being Gaussian, the other Log-Normal. In short, we have the 
following: 

itot = Il / P(Ni,Mt,i\a,k)P(M tli \Mi)dMt,i 

i 

log P(N, M\a, k) = A x log(/i) + const - a 
\ogP g (M t \Mi) = -{M t - Mi) 2 /(2a 2 M ) + const 
log P ln (M t \Mi) = -(M hl - M t j) 2 /(2(a Ml ) 2 )+ const 

(13) 

where fi = a X M + k or fi = k X M a , M t ,i stands for the M t 
for the i-th cluster and M x ,i stands for log 10 (Ma;) for the 
above mentioned form of fit. We then maximize the total 
Likelihood within a grid of resolution 100, 1000, 10000 for 
both a, k and we also vary the size of dMt,i to ensure that our 
results are robust with respect to varying grid size. We also 
try a variety of ranges for both a and k (such as < a < 20, 
-10 < a < 10 and < a < 3; and < k < 20, -20 < k < 20 
and < k < 5), and make sure we get the same answer. 
This linear fit with P g (M t \Mi) gives a = 0.455 ± 0.215 and 
k — 1.605T0.705 for 68.3% confidence intervals. This power- 
law fit with P g (M t \Mi) gives a = 0.515 ± 0.245 and k = 
1.725 ± 0.540 for 68.3% confidence intervals. See Fig. [7] for 
the data and the fit using Pi n (M t \Mi). 

3.4 N(M) distribituion: Poisson or not? 

Since we assume a Poisson distribution for N(M t ) (hereafter 
Af for simplicity in this section), we test if this is a good 
assumption by looking at ((A 2 ) - (A) 2 )/ (A). We define 
Jn = ((A 2 ) - (A) 2 )/(A), and since we only have N(M Z ), 
but not A (Number of LRG given the true measure of cluster 
mass) for each cluster, therefore, we have to consider the 
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Figure 6. The distribution of LRGs in the clusters with the fit 
to the NFW profile. Blue line: we fit the surface density of the 
LRGs (including Brightest LRGs) to a NFW profile and get a 
concentration of 17.5^ 4 'g with \ 2 = 4.29. Black line shows the 
fit of the profile of bright galaxies (normalized to this plot) from 
llLin et alj|2003) . 
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Figure 7. The number of LRG per halo N as a function of binned 
halo mass (10 14 /i^ 3 Mq), the fit (red line) is calculated by max- 
imizing the likelihood given a model of N(Mt) = a * Mt + k 
where Mt is the true measure of M200 in 10 14 /iy 3 1 MQ, assuming 
Poisson distribution of N(Mt), where Mt is the true mass of the 
cluster. We also assume a log-normal distribution for probability 
°f Pln(Mt\Mi), where Mi is the measured mass. This gives a = 
0.470 ± 0.205 and k = 1.455 ± 0.72. The blue line fit is calculated 
by maximizing the likelihood given a model of N(Mt) = k * JW™ 
and it gives a = 0.560 ± 0.250 and k = 1.575 ± 0.525. 



Figure 8. To test whether our assumption of a Poisson distri- 
bution for N(M t ) is valid, we compute {N(N - 1))/(N) 2 for the 
combined sample of 400a! + Yx. It does not deviate drastically 
from being Poisson. 



contribution of scatter from the various systematic effects 
we mentioned in i|3.1l 



(N?) (Mi ) = (N.ntf + (Nint) + Y + Z + W 
Y = 2f{N int )(N) (Mi) 

Z = J dMt[P(M t \Mi)({N) 2 (Mt) + V) 

V = -f(N) (Mt) 
W = f(l-f)(N) (Mi) (14) 

where / = f(z p ,z c ,z 3 ) and is defined in i|3.1l Ni„t is the 
number of interloper as discussed in H3.1l and X(m-) (-^(M t )) 
means the quantity X conditioned on Mi (Mt). 

Subtracting (N(Mi)) 2 M .- ) from the equation will reduce 

to: 

(N?) (Mi ) - (iV(Mi)}( Mi) =PQ + R + S + T 
P = fl+f(l-f) 

Q={N)(M i} 

R = {N int ) 



f( / dMt([P(Mt\M t )(Ny (Mt) 



T : 



-f(N) 



(Mi) 



(15) 

Note that (N) {Mi) = J dM t [P{M t \Mi){N) {Ml}) \. 

A Poisson distribution is completely characterized by its 
first moment, 7 would be 1 if the distribution is completely 
Poisson. We calculate the 7 from the combined sample of 
400d and Yx (please refer to [|4] for a description of Yx 
sample) sample and use (N) from the fit of N(M) — k x 
(M/10 14 )". We bin the cluster such that there are equal 
number of clusters in each mass bin (See Figure [8}. We find 
that 7 = 1.428 ± 0.351 and thus the N(M t ) distribution is 
consistent with being Poisson. 

Furthermore, one important ingredient of halo occupa- 
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Figure 9. To test if the satellite LRG distribution is Poisson, we 
calculate «(JV - 1) 2 > - {(N - 1)) 2 )/({N - 1)) for the combined 
sample of AOOd + Yx. It is consistent with being Poisson. 



tion distribution is the assumption of Poisson distribution of 
the satellite galaxies. We test the assumption here by com- 
puting (((N - 1) 2 > - ((N - 1)) 2 )/((N - 1)} for the (N - 1) 
distribution (see Fig. [9} in a similar way as we compute 7 for 
the (N) distribution. We find that 7^-1 = 1.823±0.496 and 
so the satellite LRG distribution in clusters is also consis- 
tent with being Poisson. However, one should note that it is 
mathematically impossible for both N and N — 1 to be both 
exactly Poisson for the same distribution. The errorbars are 
calculated via a(X) = y/[(X - X) 2 )/N c . 

3.5 LRG Multiplicity Function 

Finally we study the multiplicity of LRGs in clusters 

(Fig. EE). 

We calculate the multiplicity function by counting the 
l/Vmax weighted number of cluster in each bin. We compute 
Vmax (the comoving search volume of the cluster) by: 

i.We find the flux of the cluster via the following 
iBurenin et ail (|2006t ): 



/ 



(16) 



<±nd L (z) 2 

where L is the luminosity of the cluster, cIl{z) is is the cos- 
mological luminosity distance, K (z) is the K- correction fac- 
tor fo r X-ray clusters (for more details see IBurenin et al.l 

(Hi)). 

ii. We find the comoving search volume that each cluster 
with luminosity L can be detected given by the following: 



(17) 



where P se i(f,z) is the selection efficiency of the 400 square 
degrees ROSAT PSPC Galaxy Cluster Survey provided 
by A. Vikhlinin and R. Burenin (private communication), 
dV/dz is t he cosmologica l como ving volume per redshift in- 
terval (see IBurenin et all (|2006n for more details). 



Figure 10. The volume weighted multiplicity function of LRGs in 
these clusters. Blue (Black) line: the volume weighted multiplicity 
function for clusters with X-ray luminosities >= (<) 10 44 ergs~ 1 . 
The variance is calculated by taking the sum of l/(V^ lax ) for each 
N L ro bin. 



4 SYSTEM ATICS 

4.1 Uncertainties in the Choice of Cluster Radius 

We choose to use #200 since it is closest to the virial radius 
of the clusters. We also look at how the uncertainties of 
7*200 will affect our results. r2oo is accurate up to ~ 10% 
ijReiprich fc Bohringenl2002h . We calculate the following to 
determine the effective number of LRGs we are missing due 
to uncertainties in r2oo: 



r200 / , 1-1'"200 

p(r)dr/ I 
10 Jo 

and 



p(r)dr = 0.95 



p(r)dr/ 



0.9r 200 



p(r)dr = 1.06 



(18) 



(19) 



We set the density profile p(r) as a NFW profile with 
concentration of 8 (which is approximately what we get 
when we fit the surface density of the cluster when we ex- 
clude the BCG). This shows that the uncertainties in #200, 
thus r2oo only affect our estimation of N(M) at the level of 
~ 5%. 



4.2 Mass Estimation and Sample Selection 

Cluster mass estimation is crucial in our analysis, as it de- 
fines the cluster virial region to search for member LRGs, 
and provides a fundamental radius to scale the distance of 
LRGs to cluster center. We infer clust er mass through the X- 
ray lu minosity-mass scaling relation dReiprich fc Bohringeij 
2002) , which h as been shown as a unbiased estimator 
( Reiprichl [2006) . Compared to other X-ray-based cluster 
proxies such as temperature and Yx (the product of gas 
mass and temperature, which is proportional to the thermal 
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Figure 11. Distribution of the clusters on the mass— redshift 
plane. The squares, red points, and blue points denote the whole 
400d survey sample, the subsample used in this study, and the 
Yx sample, respectively. 



Figure 12. Combining the Yx sample (green triangles), we have 
a larger mass coverage, thus giving stronger constraints on the 
slope. We have fairly similar fits between the two different models 
(linear (red) and power-law (blue)). The fits are also consistent 
with the respective fits using only clusters from 400d survey. 



5 DISCUSSION AND SUMMARY 



energy of the cluster iKravtsov et al.l [20061 ). Lx~M correla- 
tion shows higher degree of scatter. We therefore seek for 
another cluster sample with better measured mass (despite 
without well- defined selection criter ia). 

Recently, iMaughan et al.l (|2007l ) have presented a large 
cluster sample selected from the Chandra archive, for which 
the cluster mass is inferred from Yx, and the cluster center 
is inferred from the Chandra images. 26 of these clusters lie 
within our SDSS DR5 masks and the redshift range 0.2 < 
z < 0.6. 16 of these 26 clusters do not overlap with our 400d 
sample and we use them to examine the results presented in 
§jQ£I]fc|lg] (hereafter the Y x sample). 

Because of the flux-limited nature of the 400d survey, 
low mass (~ 10 14 Mq) clusters will be only detected at lower 
redshifts. In Fig. EQ we show the mass distribution of the 
whole 400d sample (open squares) and the subsample used 
in our analysis (red points) within 0.2 < z < 0.9. It shows 
that our sample is a random subsample of the whole 400d 
sample Interestingly, at z ~ 0.3 the 400d survey clusters 
cover a larger range in mass than at other redshifts. It is 
also clear that the halo occupation numbers of clusters of 
lower masses in Fig. [7] would be biased to those at z < 0.3. 
In Fig. ED we also show the distribution of the Yx sample. 
Very curiously, the distribution of this sample on the mass- 
redshift space seems to be roughly orthogonal to that of our 
400d sample. Since our results derived from the Yx sample is 
consistent with those based on the 400d sample, we combine 
the two samples to enhance the mass coverage (especially for 
clusters at z > 0.4) and the statistical signal. We calculated 
the N(M) for the combined sample and assuming power- 
law model, we have N(M) = k * (M/10 14 ) a , where a = 
0.620 ± 0.105 and k = 1.425 ± 0.285, see Fig [H A more 
detailed result table is shown at Table 2. 



5.1 What is a good mass tracer? 

As clusters are becoming more important cosmological tools, 
we need to characterize the masses of clusters more than 
ever. "What is a good mass tracer?" has been a very well- 
motivated question. 

Here we try to investigate a few options that people 
have suggested before as possible solutions: 

Fist, as we see earlier in £13-31 the mean number of LRGs 
does not trace the masses accurately. 

We quantify this by looking at the scatter of the Nlrg — 
M relation by the following quantities in a 3 mass bins: 



a(ln(N LRG )) 



V(Nlrg) 



(20) 



a(ln(M)) = -a{ln{N LRG )) 
a 

where a is as defined in N(M) = k x (M/10 14 ) a . We found 
that the scatter in Iu(Nlrg) (ln(M)) in low, middle and 
high mass bins are 0.332 (0.535) dex (at M = 2.22 x 
10 14 h^M e ), 0.281 (0.452) dex (at M = 3.81 x lO 14 /if 3 1 M ) 
and 0.21 (0.340) dex (at M = 9.56 x lO 14 /^ 1 M ) respec- 
tively. 

Second, we look at the luminosities of the central LRG. 
As previous studies suggested in some bands, the brightest 
clust er galaxies traces the mass of the cluster l)Lin fc Mohrl 
|2004| ) and that the brightest cluster galaxies tend to be the 
central galaxies of the cluster, we look at the relation be- 
tween the luminosities of the brightest LRG in clusters and 
their X-ray masses. However, the correlation between the 
distribution of luminosities of central LRG and the masses 
of the clusters in our sample does not look promising(see 

Fig. Eg. 
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Figure 13. (Above) Luminosities of the central LRGs for each 
cluster. (Below) Luminosities of the brightest LRGs for each clus- 
ter. 



We then look at the correlation between the luminosities 
of the brightest LRG and their cluster X-ray masses. How- 
ever, it does not se em to be promising either (see Fig. [T3)) . 
This is also seen in iLin fc Mohrl l|2004 ) when we look at the 
same mass range and when one looks at the correlation be- 
tween t he brightest LRGs an d the richness of the maxBCG 
catalog |Koester et al.ir~2007l ) . there is not a strong correla- 
tion for 14,000 clusters (R. Reyes 2007, private communica- 
tion). However, there are several caveats that would require 
further investigations, such as the possibility of photo-z fail- 
ure for the CLRG or BLRG in the clusters and possible 
photometry problem that could destroy the correlation. We 
look into the available spectroscopic data in SDSS and found 
no extra LRGs that are targeted by the SDSS spectroscopy. 
This rules out the possible missing LRGs that have M r of 
range ~ -20.8 (at z = 0.2) and ~ -22.5 (at z = 0.6). Fur- 
thermore, as we investigated earlier, only 4 clusters do not 
have LRGs and we find ~ 70% of BLRGs lie in the central 
~ 20% of the virial radius, thus, most clusters do have a 
LRG at their centers. If we are missing Brightest LRGs in 
centers of clusters, we need to expect the scenario of hav- 
ing more than 1 LRG at the central ~ 20% of cluster virial 
radius to be prevalent. This scenario is not supported by 
the distribution of LRGs as shown in Figure [4] Given the 
caveats and findings here, we conclude that further work 
will be needed to make this more quantitative, especially to 
quantify the effect of photometry errors on the correlation. 



5.2 Evolution of Massive Galaxies 

Insights into the evolution of massive galaxies in clusters 
may be gained by comparing some of the results presented 
in Sj2] with the properties of cluster galaxies in the local 
Universe. For consistency with our LRG selection, we se- 
lect nearby cluster galaxies by the requirement that they 



are one magnitude more luminous than the characteristic 
magnitude (M* — 1). 

We first examine t he spatial dis t ribut ion of massive 
galaxies within clusters. ILin fc Mohrl (I2004T ). with a large 
sample of clusters at z < 0.2, find that luminous cluster 
galaxies (Mk < —25) follow an NFW profile with concen- 
tration of 18.2 (5.8), when the brightest cluster galaxy is 
included (excluded). This result is in very good agreement 
with our finding in H3.2I 

The second comparison is made with the halo occu- 
pation number. We construct the occupation number for 
Mk < —25 wit h the z < 0.1 cluster sample presented in 
ILin et all (|2006l ). No color selection analogous to that pre- 
sented in £|2.2I is used. However, for such a luminous mag- 
nitude range, the contamination of blue galaxies should be 
minimal. Nevertheless, we should regard the mean occupa- 
tion number N thus obtained as a upper limit. We find that 
the nearby N—M relation is similar to that shown in ij3 .3 1 

Taken at face value, these comparisons seem to sug- 
gest that there is not much evolution in the massive galaxy 
populations between z ~ 0.5 and z w 0. The occupa- 
tion number comparison basically suggests that the shape 
of the luminosity function is similar in clusters at these 
two epochs, after the passive evolution has been taken into 
account. This is consisten t with several previous studies , 
both for cluster galax i es (ILin et al.l 20061; lAndreonl 120061 ; 
|Pe Propris et al.l 120071; iMuzzin et al.l |2007|) and the field 



population e.g. dWake et al.ll2006l ; iBrown et alll2007r i. 

However, evidence for mergers that produce massive 
galaxies has been found (|van Dokk um 2005; Bell et al.ll2006l ; 
some other references). In the ACDM model, formation of 
massive objects through mergers of less massive ones is a 
generic feature. To reconcile the apparent no-evolution of the 
aforementioned bulk properties with this picture, we suggest 
two considerations. (1) Irrespective of the role of mergers in 
the formation and evolution of the LRGs, their spatial dis- 
tribution seems to be similar ou t to z ~ 0.5 . This is similar 
to the "attractor" hypothesis of lGao et all (|2004l ). (2) The 
degree of evolution, be it an increase in the number of LRGs 
due to mergers of the host halo with less massive halos, or 
a decrease due to dynamical processes (e.g. tidal disruption, 
mergers) , would be seen more clearly through (Monte Carlo) 
simulations where the merger history of the halos is fully fol- 
lowed. In a companion paper s uch an approach is adopted t o 
infer the merger rate of LRGs IjConrov. Ho fc White 1 12007). 

5.3 Summary 

We investigate statistical properties of LRGs in a sample 
of X-ray selected galaxy clusters at intermediate redshift 
(0.2 < z < 0.6). The LRGs are selected based on carefully 
designed color criteria, and the cluster membership is as- 
sessed via photometric redshift. We put constraints on spa- 
tial distributions of LRG within clusters, namely the radial 
distribution. We find that the distribution of brightest LRGs 
in c luster to be concentrated as discussed in previous stud- 
ies (I Jones fc Formadll984l ; ILin fc Mohrll2004D . We also find 
that the radial distribution can be fit by a NFW profile with 
a concentration of 17.5l 4 'g with \ 2 = 4.29 when we include 
the brightest LRG. When we do not include the brightest 
LRG, we find concentration of 6.0l j'9 with \ 2 = 6.6. Consid- 
ering the sample size and mass errors on our sample, we use 
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the maximum likelihood method to find the best fit param- 
eters for halo occupation distribution (N(M)). The result 
depends on what kind of models we adopt, but are fairly 
insensitive to what model we use, results are shown in Ta- 
ble 2. 

Uncertainties in photometric redshifts are taken into 
account by including different possible effects such as in- 
terlopers and missing LRGs due to errors in photometric 
redshifts (see SJ3]). We estimate that the errors in cluster ra- 
dius can only contribute to our uncertatinty in N(M) at 
the level of ~ 5%. Errors in mass estimation are fully taken 
into account throughout the analysis. We also employ an in- 
dependent sample of better measured masses (Yx sample) 
to test the mass estimation of our sample. However, we do 
implicitly assume that the scatter of Al — Lx relation does 
not correlate with N(M) during the analysis. The result we 
derive from a combined analysis of both sample on N(M) 
is consistent with using our sample alone (see Table 2). We 
also find that there are no obvious good mass tracer as we 
look at different correlations between various quantities of 
clusters and their galaxies. Last, we discuss the evolution 
of massive galaxies from different perspectives. We conclude 
that it would be important to study low-z LRG p opulation 
to bet ter constrain the evolution of the population (|Ho et al.l 
I 20071 ). 
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Table 1. Basic parameters of our cluster sample. Naming scheme 
follows the cluster number as given in ta ble 4 o f Burcn n et al.l 
ll2006ll those start with s are from table 5 of lBurenin et al. N2006I) . 
and are not part of the main sample of serendipitous 400d survey. 
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N(M) 




Data 




Poisson+Gaussian 


Poisson+log-Normal 


a * M 4 


k 


400d 




a = 0.455 ± 0.215 


a = 0.470 ±0.205 
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400dH 


-Y x 


a = 0.460 ± 0.090 


a = 0.430 ± 0.085 










fc = 1.500 ±0.435 


k = 1.515 ±0.435 


k*M a 




400d 




a = 0.515 ± 0.245 


a = 0.560 ±0.250 










k = 1.725 ±0.540 


k = 1.575 ±0.525 


k * M a 




400dH 


-Y x 


a = 0.625 ± 0.110 


a = 0.620 ±0.105 










k = 1.470 ±0.285 


k = 1.425 ±0.285 



Table 2. Results of maximizing likelihood assuming different pa- 
rameters. This table describes the model of the N(M) in the first 
column, dataset we use in the "Data" column, results of maximiz- 
ing likelihood by assuming Poisson distribution for N(M) (both 
column 3 and 4), Gaussian and Log-Normal distribution for the 
cluster mass distribution for column 3 and 4 respectively. 
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